Indication of a Co-Existing Phase of Quarks and Hadrons in Nucleus - Nucleus 

Collisions 
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The variation of average transverse mass of identified hadrons with charge multiplicity have been 
studied for AGS, SPS and RHIC energies. The observation of a plateau in the average transverse 
mass for multiplicities corresponding to SPS energies is attributed to the formation of a co-existence 
phase of quark gluon plasma and hadrons. A subsequent rise for RHIC energies may indicate a 
deconfined phase in the initial state. Several possibilities which can affect the average transverse 
mass are discussed. Constraints on the initial temperature and thermalization time have been put 
from the various experimental data available at SPS energies. 
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Results based on QCD (Quantum Chromodynam- 
ics) renormalization group approach predict that 
strongly interacting systems at high temperatures 
and/or densities are composed of weakly interacting 
quarks and gluons yj, [2j due to asymptotic freedom 
and the Debye screening of color charge. Nucleus- 
Nucleus (A- A) collisions at very high energies may 
create situations conducive for the formation of a 
thermodynamic state where the properties of the 
system are governed by quarks and gluonic degrees 
of freedom. Such a state is called quark gluon plasma 
(QGP). QCD lattice gauge theory suggest that the 
critical temperature (T c ) for such a transition is 
- 170 MeV 0. This has led to intense theoret- 
ical and experimental activities in this field of re- 
search 0. 

Experimental detection of the QGP in A-A colli- 
sions is a non-trivial task because of the small space- 
time volume of the system. Among many signals of 
the QGP formation 5], one of the earliest is based 
on the relation of the thermodynamical variables, 
temperature and entropy to the average transverse 
momentum and multiplicity, respectively. This was 
originally proposed by Van Hove in the context of 
proton-proton collisions @- It was argued that a 
plateau in the transverse momentum beyond a cer- 
tain value of multiplicity will indicate the onset of 
the formation of mixed phase of QGP and hadrons; 
analogous to the plateau observed in the variation of 
temperature with entropy in a first order phase tran- 
sition scenario. One hence looks for the variation of 
average transverse momentum ((pt)) or transverse 
mass ((tot), tot = \fpr + to 2 ) with respect to the 
total number of particles produced per unit rapid- 
ity (dN/dY) in A-A collisions at high energies. An 
increase in energy density should increase (%) till 
it reaches a critical density for phase transition to 
occur, where it should then show a plateau due to 
formation of mixed phase. Further increase of energy 
density should again increase the (tot) of produced 
particles. 



Several attempts have been made so far to look 
for such signals. While the data from cosmic ray 
experiments are inconclusive, the data from A-A 
collisions in the laboratory for fixed center of mass 
energies (y/s) with different centralities (impact pa- 
rameters) H U E3 have so far not shown this kind 
of behaviour. This implies that a mere change in 
the centrality of the collisions does not change the 
energy density of the system formed after the col- 
lisions required to create any change of phase. So 
it is imperative to study the variation of (tot) with 
dN c h/dY for a broad range of beam energies for fixed 
centrality. 

In this letter the variation of (tot) with dN c h/dY 
is examined for AGS, SPS and RHIC energies span- 
ning s/s from 2A GeV to 200A GeV. We then carry 
out explicit theoretical calculations to understand 
the observed behaviour in terms of the properties of 
the matter produced in the initial stages. 

In fig. H the variation of (tot) with charge mul- 
tiplicity is depicted for pions, kaons and protons at 
AGS, SPS and RHIC energies around mid-rapidity. 
The data shown here correspond to central events 
for different colliding systems, (produced) particle 
types and center of mass energies (see Table D). The 
experimental data on transverse mass spectra can be 
parametrized as, 
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where the inverse slope parameter T e f f is the effec- 
tive temperature (effective because it includes the 
contribution from both the thermal and collective 
motion in the transverse direction). The variation 
of T e ff with y/s for kaons has recently been shown 
in Ref [llj . The average transverse mass of the par- 
ticles obtained from Eq. ^is: 
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FIG. 1: Variation of (mr) with produced charged parti- 
cles per unit rapidity at mid rapidity for central collisions 
corresponding to different y/s spanning from from AGS 
to RHIC. The error bars reflect both the systematic and 
statistical errors in obtaining T e //. 



TABLE I: The details of the experimental data used 
the analysis along with the references. 
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From the results shown in ng^ one observes an 
increase in (m^) with dN c h/dy for AGS energies fol- 
lowed by a plateau for charge multiplicities corre- 
sponding to SPS energies for all the particle types, 
pions, kaons and protons. This may hint at the pos- 
sible co-existence of the quark and hadron phases. 
For charge multiplicities corresponding to RHIC en- 
ergies, the {tut) shows an increasing trend indicat- 
ing the possibility of a pure QGP formation. 

It is essential now to investigate whether the above 
experimental variation of (my) can arise from vari- 
ous physical effects other than due to QGP forma- 
tion. We will address this issue by trying to answer 
the following questions : (a) How does (3+1) dimen- 
sional hydrodynamical evolution of the system, with 
and without QCD phase transition, formed in heavy 
ion collisions affect the {tut) at freeze-out ? (b) Does 
the analysis of experimental data on hadrons, pho- 
tons and dileptons indicate similar values of the ini- 
tial temperature (T^)? How does the value of Tj com- 
pare with the T c predicted by lattice QCD Q? (c) 
How does the (tot) get affected by the gain in trans- 
verse momentum by hadrons through their succes- 
sive collisions (Cronin effect [2fJ) with the particles 
in the system ? (d) How does the (m<r) of produced 
particles vary with yfs in a standard event generator 
based on the principle that nucleus-nucleus collision 
is a superposition of nuclcon-nuclcon collisions ? 

(a) To understand the variation of (my) with 
dN c h/dY for a system undergoing a (3+1) di- 



mensional expansion we solve the hydrodynami- 
cal equations for various initial energy densities 
( e i) uM- The initial radial velocity is taken as zero. 
The momentum distribution of particles originating 
from a system undergoing transverse expansion |2l| 
with boost invariance along the longitudinal direc- 
tion [22, contains the effect of thermal motion as 
well as transverse flow. Hence it is imperative to 
examine how the value of the inverse slope or aver- 
age transverse mass is affected by transverse motion. 
To achieve this we proceed as follows. We evaluate 
the mr spectra for various hadrons at the freeze-out 
point as a function of for two different kind of 
space time evolution scenarios. 

Firstly, we consider a hadronic gas at various val- 
ues of €i from ~ 0.05 to 4 GeV/fm 3 without intro- 
ducing a QCD phase transition. For all these ini- 
tial energy densities the (3+1) dimensional hydrody- 
namic equations are solved to obtain the freeze-out 
surface at a temperature of 120 MeV. This is used as 
an input to evaluate the transverse mass spectra of 
hadrons. The hadronic equation of state (EOS) |2^] 
contains hadronic degrees of freedom up to strange 
sector 24]. From the nix spectra of hadrons (pi- 
ons, kaons and protons) corresponding to each q we 
evaluate the (m^)- This is plotted in Fig. |2 (dashed 
line) . 

In the second scenario, we introduce the QGP 
(composed of up, down, strange quarks and glu- 
ons) at an energy density above £q ~ 1.7 GeV/fm 3 
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FIG. 2: Variation of (mr) with initial energy density 
(ei) for a QGP scenario (solid line) and hadronic scenario 
(dashed line) obtained from (3+1) dimensional frydrody- 



corresponding to critical temperature, T c = 170 
MeV which converts to the hadronic phase at an 
energy density, eh ~ 0.56 GeV/fm 3 at the same 
temperature. The effective degeneracy (g e ff) at T c 
for the hadronic phase is ~ 16 obtained from the 
study of hadronic EOS in |2^|. So in this sce- 
nario we have a mixed phase within the energy range 
0.56 < e(GeV/fm 3 )< 1.7. For e > 1.7GcV/fm 3 a 
QGP bag model EOS has been used. The variation 
of rriT as a function of for this scenario is shown 
in Fig. [21 (solid line). It is clear that (tot) varies at a 
much slower rate for the second scenario. Although 
a complete plateau structure is not observed, it can 
be argued that the formation of a mixed phase slows 
down the rate of increase of (mr) with multiplicity. 
In fact, the slope of the curve in the case of pions 
for the first scenario is larger than the second by a 
factor of 2.5 at e, = 1 GeV//m 3 . The (m T ) for pi- 
ons is about 400 MeV for values of e< corresponding 
to mixed phase, which is similar to those obtained 
experimentally. 

(b) The lower bound on the T$ realized at SPS 
energies can be obtained from the measured (tot) 
provided we can estimate the effects of the flow. The 
inverse slope of the tot spectra contains the effect 
of flow and can be parametrized as |25|. 
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where (v r ) is the average velocity and T — Tp is the 
"true" freeze-out temperature. We start by noting 
that (tot) is actually the average energy, (E) of the 
particle at zero rapidity (E = mp coshy) and in a 
thermal system the average energy is proportional to 
temperature. Therefore, for a system which has zero 
transverse velocity but non-zero longitudinal veloc- 
ity one must have 



(to t ), > (tot, 



(4) 



because of the loss of thermal energy due to work 
done in longitudinal expansion. Here the subscript 
i (f) refers to initial (freeze-out) state. Assuming 



TO T/ 



mr} expt , (i.e. all the hadrons are emitted 



from the freeze-out surface) and using Eq's.[21and2| 

we get 
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A value of (v r ) ~ 0.4 reproduces the slope of the tot 
distribution of pions at SPS energies reasonably well. 
For this value of (v r ) we get Tp ~ 120 MeV from 



Eq.|3 indicating, (mr) 



, + T F + T F /(m + T F ) 



315 MeV. It should be noted here that the effects of 
the transverse flow is "subtracted" out from (tot) 
now. Using Eq.|3]we obtain T t > 158 MeV. 

Photons and dileptons emitted from the strongly 
interacting matter constitute one of the most effi- 
cient probes to measure the initial temperature of 
the system It was shown earlier [2(| that di- 
rect photons (measured by WA98 collaboration |27jl 
having transverse momentum in the range 1 < pp 
(GeV)< 2.5 is expected to originate from a thermal 
source. The inverse slope of this spectra is ~ 275 
MeV. For (v r ) ~ 0.4 we obtain an average temper- 
ature T av ~ 178 MeV by using Eq. [3J A similar- 
value of T av is obtained from the pp distribution of 
e + e~ [28| . Since photons and dileptons are emitted 
from all the stages of space time evolution, T av satis- 
fies the condition, Ti > T av > Tp. All these indicate 
that a value Ti > 170 MeV may have been realized 
at SPS energies. 

The total multiplicity of hadrons is related to the 
initial temperature (Ti) and thermalization time (t$) 
as: 



dN 
dY 



^9eff^R 2 T t 3 r t /3.6, 
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R is the radius of the colliding nuclei and g e f / is the 
effective statistical degeneracy. Taking dN/dY ~ 
700 corresponding to the highest SPS energies and 
using Eq. 0we get n < 1.2 fm/c ~ 1/Aqcd where 
Aqcd is the QCD scale parameter. 

(c) Recently, there have been some attempts to ex- 
plain the transverse mass spectra of produced par- 
ticles without introducing a transverse expansion. 
The hadrons gain transverse momentum through 
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successive collisions with the particles in the sys- 
tem. This can be treated as a random walk of the 
hadrons p9| in the medium. The Tot spectra of pi- 
ons at SPS energies can be reproduced by this model, 
but it fails to reproduce the slope of kaon and proton 
distribution at the same colliding energies. 

(d) As a last consideration, we investigate how 
the (tut) of produced particles vary with y/s of A-A 
systems in a standard event generator based on the 
principle that the A-A collision is a superposition 
of nucleon-nucleon collisions. This is done by using 
the event generator, HIJING 30] which contains the 
effect of mini-jets. The results of this calculation do 
not show the variation of (tot) with ^fs as has been 
observed experimentally. 

In summary, a detailed analysis of the experimen- 
tal data of the transverse mass spectra of identi- 
fied hadrons, at AGS, SPS and RHIC energies, show 
a characteristic behaviour of (mj) as a function of 
multiplicity similar to that proposed by Van Hove. 
The plateau in (tot) corresponding to the multiplic- 
ity at SPS energies may hint at the formation of 
a mixed phase of QGP and hadrons in the initial 



stages. The average temperature obtained from the 
analysis of photons and dileptons lends support to 
this. We have shown by solving the hydrodynam- 
ical equations that the presence of a mixed phase 
really slows down the growth of (tot) with dN/dY 
substantially. In the present work all the hadronic 
degrees of freedom up to mass 2.5 GeV has been in- 
cluded in the EOS; a fewer hadronic degrees of free- 
dom will slow down this growth even more. Since 
a lower bound of the initial temperature is obtained 
here, the formation of QGP at a higher temperature 
cannot be ruled out. However, even if QGP phase is 
created at SPS, its space time volume could be too 
small to affect the dynamics of the system and hence 
our conclusions. The upward trend in the slope of 
(tot) at RHIC energies is indicative of a deconfined 
phase in the initial state. 
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